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ABSTRACT 

We measure the cosmological matter density by observing the positions of baryon acoustic oscil- 
lations in the clustering of galaxies in the Sloan Digital Sky Survey (SDSS). We jointly analyse the 
main galaxies and LRGs in the SDSS DR5 sample, using over half a million galaxies in total. The 
oscillations are detected with 99.74% confidence (3.0c assuming Gaussianity) compared to a smooth 
power spectrum. When combined with the observed scale of the peaks within the CMB, we find a 
best-fit value of £Im = 0.2561qq24 (68% confidence interval), for a flat A cosmology when marginal- 
ising over the Hubble parameter and the baryon density. This value of the matter density is derived 
from the locations of the baryon oscillations in the galaxy power spectrum and in the CMB, and 
does not include any information from the overall shape of the power spectra. This is an extremely 
clean cosmological measurement as the physics of the baryon acoustic oscillation production is well 
understood, and the positions of the oscillations are expected to be independent of systematics such 
as galaxy bias. 

Subject headings: cosmology: cosmological parameters, large-scale structure of universe 



1. INTRODUCTION 



Baryonic Acoustic Oscillations (BAO) are predicted 
in the matter distri bution wit h a calibration that 
depends on n M h 2 (Isiikl H96l iPeebles fe Yi 
Sunvaev fc Zel 'dovich 19701 : iBond fc Efstathia 
1987; Holt zmanl 119891 ). The oscillations arise because 
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sound waves in the coupled baryon-photon plasma af- 
ter an inflationary epoch will lead to the expansion 
of the baryonic material in a spherical shell around a 
small perturbation, reaching a radius rs(z*), the co- 
moving sound horizon size at recombination, before 
sound waves are no longer supported within the plasma 
(jBashinskv fc Bertschingeill200lLl2002h . At the high red- 
shifts of interest the vacuum energy can be neglected, and 
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rs( z *) can be simply written (|Hu fc Sugivamal[l995f ) 

1 



r s (z,)/fe Mpc = 



100f^{ 2 Jo (« + acq) 17 

The expansion factor a = (l + z) -1 and a 
values at recombination and matter-radiation equality 
respectively. Thus rs(z*) depends on the matter den- 
sity fijvf through the expansion rate and the recombina- 
tion redshift. For a baryon density ^h 2 ~ 0.02, we can 
approximate cs — 0.90C/V3- Inserting = 1100 and 
a cq = (23900 VL M h 2 )- 1 gives r<?(z«) = 109/i- x Mpc for 
VL M = 0.24 and h = 0.73. 

In real space this leads to a peak in the correlation 
function at rs(z*). In Fourier space, this process leads to 
oscillations in the power spectrum in the same way that 
the transform of a top-hat function yields a sine function. 
The wavelength of these oscillations for Qm — 0.24 and 
h = 0.73 is k s = 2tt/109 = 0.06/iMpc" 1 . Numerical 
simulations have shown that a number of subtle correc- 
tions are required to this simple picture, although these 
corrections do not significantly affect the underlying im- 
portant physics. On large scales, there is a phase shift in 
the position of the oscillations d ue to the contribution o f 
the baryons to the drag epoch (Eisenste in fc Hulll998f) . 
We must also consider the d amping of the oscillations on 
small scales at high redshift (Silk 1968), and, when mod- 
elling the oscillations at low reds hifts, due to structure 
formation ijEisenstein et al.ll2006f ). 

The resulting acoustic peaks in the cosmic microwave 
background (CMB ) have now been o bserved with 
extre me precision ()Hinshaw et al.l 120061 : ISnergel et al.1 
l2006h . However, interpreting the observed angular sepa- 
ration of these peaks in terms of the physics of the early 
universe requires knowledge of the angular diameter dis- 
tance to the last scattering surface - in the context of flat 
A cosmological models this leads to a parameter degener- 
acy between the matter density Qm and the Hubble con- 
stant h, and models with the same value of Vt°^ 275 h have 
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the s ame projecte d acoustic horizon s cale (IPercival et all 
120021: iPage et all 120031 : ISpergel et al.ll2006h (throughout 
this paper, a "flat A" cosmological model implies a 
universe that is spatially flat with a time-independent 
"dark energy" component, i.e., a cosmological constant). 
Theoretically, we should expect th e oscillations to sur- 
vive in the galaxy power spectrum dMeiksin et aL 
Springel et altooa ISeo fc Eisensteinll2005l ; IWhit" 
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Springel et alJl 2005: Sc o fc Eisensteinll2005l : IWhitdl2005i 
Eisenstein et al.ll2006f ). and the effects of baryons have 
been previously detected on large scales in the cluster- 
ing of galaxies (IPercival et alj [200lt iMiller et all 120011 : 
ICole et all 120051: lEisenstein et all 120051: iHuetsil l2006f ). 
These low redshift observations have a different depen- 
dence on the cosmological distancc-redshift relation com- 
pared with the CMB because of the different angular and 
radial projection of these features. By comparing the two 
observations we can therefore probe the cosmological ex- 
pansion history in addition to the physics of the BAO 
production. In this paper we use the BAO measured in 
the SDSS DR5 galaxy sample to set tight, clean con- 
straints on the cosmological matter density. 

2. MEASURING THE SDSS POWER SPECTRUM 

The data and method used for calculating the 
redshift-space power spectrum o f the latest SP S S sam - 
ple are described in detail in IPercival et al.1 (2006). 
In this section we summarise this information pay- 
ing particular attention to impo rtant detai ls for 
the BAO detection. The SPS S (lYork et all 
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lease 5 (PR5) sample represents the largest volume of 
the Universe that has been mapped to date. In total we 
have 522280 galaxy spec tra, with 465789 of those spec- 
tra being main galaxies (jStrauss et al.ll200"2l ) selected to 
a limiting magnitude r < 17.77, or r < 17.5 in a small 
subset of the early data from the survey. The remain- 
ing 56491 galaxies ar e Luminous Red Galaxies (LRGs; 
lEisenstein et al.l 1200 ll ). which form an extension to the 
survey to higher redshifts 0.3 < z < 0.5; this extension 
covers most of the volume mapped. Of the main galaxies, 
21310 are also classified as LRGs, so our sample includes 
77801 LRGs in tota l. In this paper, and its companion 
(|Percival et al.ll2b06l ). we analyse the combined sample of 
main galaxies and LRGs thereby including correlations 
between the two samples. Although the main galaxy 
sample contains significantly more galaxies than the LRG 
sample, the LRG sample covers more volume, and there- 
fore contains almost all of the cos mological information . 
However, as shown in figure 8 of IPercival et~aT1 (|2006h . 
the contribution of pairs of galaxies where one galaxy is 
a main sample galaxies and one is a LRG are not negligi- 
ble compared with LRG-LRG pairs, justifying the added 
complexity of an analysis using both data sets. 

The sample now contains 60% more LRGs than con- 
sidered in th e first measurement o f BAO in the SPSS 
LRG sample (Eise nstein et al.l [20051 ): with the increased 
precision due to this increase in volume and the main - 
LRG pairs of galaxies, we can now attempt to derive 
constraints on the matter density that do not rely on 
the overall shape of the Cold Park Matter (CPM) power 
spectrum but only on the peak locations relative to a 



smooth underlying spectrum. Although ignoring the 
overall shape removes information, we gain in robust- 
ness due to the increased simplicity of the physics pro- 
ducing the cosmological constraint; measurements of the 
matter density from the overall shape of the galaxy 
power spectrum are potentially affected by galaxy bias 
(|Percival et al.ll2~006l ). or the form of the fluctuation spec- 
trum from the inflationary model. 

The redshift-space clustering power spectrum of this 
sample has been calculated using a Fourier base d tech- 
nique ((Mdmi^iiiniilll Hird^ririn[200l). This 
method uses a simple model for the relative bias of galax- 
ies to remove effects in the power spectrum due to pairs 
of galaxies with different expected clustering amplitudes. 
The power spectrum is recalculated for 31 flat A cosmo- 
logical models, with matter density 0.1 < VLm < 0.4 and 
AH,m = 0.01. For each of these cosmological models, 
we have created 200 log-normal cat alogues (using the 
method described in lCole et al.ll2005l ) with power spec- 
tra calculated using a linear CPM model with parameters 
chosen to approximately match the amplitude and shape 
of the recovered power for 0.01 < k < 0.15 /iMpc -1 . 
These power spectra are used to calculate a covariance 
matrix for each model, although we then minimise the 
effect of Monte-Carlo noise in each element in the set 
of covariance matrices as a function of £Im by smooth- 
ing using a 4-nodc cubic spline with nodes at flm — 
0.1,0.2,0.3,0.4. This ensures a smooth progression in 
the error estimation over the set of cosmological mod- 
els. The convolving effect of the survey geometry on the 
power spectrum has also been quantified by spline fitting 
the Fourier transform of the window function. This fit 
is used to smooth all models before comparison with the 
data. 

3. MODELLING THE BARYON ACOUSTIC OSCILLATIONS 

We model the true galaxy power spectrum on scales 
< k < 0.3ft.Mpc _1 with a two component model. 
The overall shape was matched using a cubic spline 
fit with 8 nodes separated by AA; = 0.05 /iMpc -1 and 
0.025 < k < 0.375 /iMpc -1 , and an additional node at 
k = 0.001 /iMpc" 1 . This smooth model was then mod- 
ulated by a higher frequency component, constructed as 
follows: the sinusoidal BAO term in a standard CPM 
transfer function was estimated for the parameters to 
be tested from nume rical fits (Eis enstein fc Hul 119981 ; 
lEisenstein et al.ll2006l ). including a damping term to ap- 
proximately correct for non-linear structure formation, 
and the multiplicative effect on a CPM power spectrum 
was isolated. This multiplicative term was then ap- 
plied to the smooth cubic spline power spectrum rather 
than a CPM model. The model power spectra are ad- 
justed for the effects of the survey geometry by convolv- 
ing with the appropriate window function and correcting 
for our lack of knowledge about the true mean density of 
galaxies by subtracting a multiple of the Fourier window 
function from the m odel power so that P(0) = (see 
IPercival et aT1l2006l for details). 

This procedure separates the physics of the BAO from 
that governing the overall shape of the power spectrum, 
including both cosmological and galaxy formation effects, 
and ensures that the cosmological constraints presented 
in this paper only come from the BAO and not from the 
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additional physics encoded in the power spectrum. In 
particular, forming a model power spectrum in this way 
allows for non-linear effects and galaxy bias to change the 
overall shape of the power spectrum, and damp the os- 
cillations on small scales. The model does not allow the 
BAO scale to change, although it is worth emphasising 
that the model can lead to an apparent change in the po- 
sitions of the peaks and troughs in the model power spec- 
tra, caused by the addition of a tilted smooth power spec- 
trum component (for example the 1-halo term in the halo 
model), or multiplication by such a smooth component. 
Recent theoretical models of BAO in galaxy power spec- 
tra find such an ap parent change in the po s itions of the 
peak s and troughs (|Seo fc Eisensteinl 120051 : ISmith et al.1 
l2006f) . but have not shown evidence for a change in the 
BAO scale. Because we split into the BAO and a smooth 
component, any such apparent shift in the observed BAO 
scale is removed by our analysis method. Extra multi- 
plicative or additive low frequency power spectrum com- 
ponents will change the BAO damping, which will be 
a function of the luminosity of the galaxy sample. For 
a combined sample of galaxies with different luminosity 
such as analysed in this paper, the form and amplitude 
of the small-scale BAO damping will depend on the de- 
tails of the galaxy sample. However, the damping does 
not have a significant effect on our results, so our results 
are expected to be robust to such complexities. 

Because of the large volume observed, the Fourier 
modes are not strongly correlated - correlations between 
modes drop to < 0.33 for Ak > 0.01 /iMpc -1 . The cor- 
relation scale is therefore significantly smaller than the 
BAO wavelength. However, we calculate the likelihood 
of the data given each model to be tested assuming that 
the data are drawn from a multi-variate Gaussian dis- 
tribution with covariance matrix calculated as described 
in Section [5] (including the effect of the change in the 
determinant as a function of fijvf). For each cosmolog- 
ical model tested, we calculate the maximum likelihood 
obtained after varying the values of the smooth spline 
function at the nine nodes for models where we include or 
exclude the BAO model. Ideally, excluding the BAO, we 
would find that the likelihood does not change with the 
cosmological model. However, because the data change 
with Qm we find small variations in the likelihood even 
when only fitting a smooth curve to the power spectra. 
To remove this "noise" , we only consider the ratio of the 
different likelihoods including and excluding the BAO 
signal. 

This method relies on the spline curve being able to fit 
the power spectrum shape. In order to test this, we have 
fitted two sets of 100 l inear CDM power spe ctra (cal- 
culated using the fits of lEisenstein k. Hulll998l assuming 
h = 0.73, fijif = 0.24 and n s = 0.96) with our combined 
spline+B AO model. The first set of mock data contained 
BAO with flb/flM = 0.174, while the second set had no 
BAO (fib = 0). Noise was added to the mock power 
spectra drawn from a multi-variate Gaussian distribu- 
tion with covariance matrix matched to the real data, 
and the power spectra were convolved with the window 
function of the survey (for a cosmological distance model 
with f2j\/ = 0.24). For the combined splinc+BAO model 
that we fit to these data, we have assumed that there 
is no small scale damping, and have fixed the cosmolog- 
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Fig. 1. — Results from fitting to 200 mock linear CDM power 
spectra calculated assuming h = 0.73, tljf = 0.24 and n s = 0.96, 
with our two component spline + BAO model. We plot the average 
ratio between input power and spline fit from f 00 mock CDM power 
spectra with no baryon oscillations (open circles with la errors), 
and 100 mocks with baryon oscillations assuming Q^/Qm = 0.174 
(solid circles with la errors). The expected residual in each case is 
shown by the solid lines. From these fits, we find that the differ- 
ence between input and recovered power AP(fc)/P(fc) < 0.015 for 
0.01 < k < 0.3 h Mpc" 1 , a level well below the current experimen- 
tal error (grey shaded region). 

ical parameters at the input values. The spline curves 
are then allowed to vary to fit the data, and the average 
residuals between the mock power and the spline fit are 
plotted in Fig. [TJ compared to the expected residuals. 
The results match those expected at a level well below 
the error on the power spectrum; this, and the accept- 
able average % 2 values of the fits (64.2 and 66.2 for the 
sets of models with and without BAO, given 64 degrees- 
of-freedom) show that the spline+BAO model can match 
the features expected in a linear CDM power spectrum. 

Finally, we end this section with a brief discussion on 
the node separation chosen for the spline fit. If too many 
nodes are chosen, then the spline fit can itself match the 
BAO, leading to a small likelihood ratio between models 
with and without BAO even if BAO are strong in the 
data. Conversely, if too few nodes are chosen, the ad- 
dition of BAO to the model will not necessarily match 
high-frequency features in the observed power spectrum, 
and might instead simply help to fit the overall shape. 
The node spacing adopted in this work was carefully cho- 
sen based on the analysis of mock power spectra, and 
from the results of the fits to the SDSS power spectra 
- with the chosen separation no evidence was observed 
for low frequency residuals from any of the SDSS power 
spectra fits, with or without BAO. 

4. RESULTS 

In Fig. [2 we plot measured power spectra determined 
assuming 5 different values of the matter density, divided 
by the best-fit smooth cubic spline fit (solid circles). The 
spline fits were calculated after fitting the data including 
a possible BAO contribution with fixed fib/ 'fin ~ 0.17 
and h = 0.73, and the appropriate SI a/- These data 
are compared with the BAO model, and show that the 
model and data match only if 0.2 < flM ^ 0.3, assuming 
O&AW ~ 0.17 and h = 0.73. For fl M = 0.26, the baryon 
oscillations are required with a likelihood ratio 2 In C = 
9.02, corresponding to a 99.74% confidence of detection 
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Fig. 2. — The ratio of the power spectra calculated from the 
SDSS to the smooth cubic spline fit that we use to model the over- 
all shape of the measured power spectra (solid circles with 1-a 
errors). Data are plotted using five flat A cosmological models to 
convert from redshift to comoving distance, with matter densities 
given in each panel. For comparison, in each panel we also plot the 
BAO predicted by a CDM model with the same matter density, 
h = 0.73, and a 17% baryon fraction (solid lines). The dashed 
lines show the same models without the low-redshift small-scale 
damping term. As can be seen, the observed oscillations approxi- 
mately match those predicted by this model for 0.2 < Ojvf < 0.3, 
but fail for higher or lower matter densities. 

(3. Oct assuming Gaussianity). 

Because the amplitude of the BAO depends on the 
fraction of matter that is baryonic, the SDSS observa- 
tions constrain the baryon fraction, although this con- 
straint is weak compared to constraints on £!t,/i 2 from 
CMB observations. The scale of the BAO depends on 
Qjf and h, so we need one other piece of information 
about a combination of these parameters to break this 
degeneracy and measure Qm- We consider three op- 
tions, yielding the three sets of likelihood contours in 
the h — fijvf plane that are shown in Fig. [3] First, we 
combine our low redshift BAO measur ement with the 
HST Key Project (jFreedman et al.ll2001h estimate of the 
Hubble parameter h = 0.72 ± 0.08 (Ict errors). This 
gives £Im = 0-256^q'q29 ( a ^ the error bars quoted in this 
paper span the 68% confidence interval). Here we have 
marginalized over the uncertainty in h (assumed to be 
Gaussian) and marginalised over a uniform prior on the 
baryon density 0.008 < f^/i 2 < 0.034. This estimate of 
Am assumes a flat A cosmological model, but dark en- 
ergy and spatial curvature affect the result only through 
their influence on the distance-redshift relation for the 
galaxies within our sample. Measurements of Type la su- 
pernovae already demonstrate that this relation is close 
to that expected for a flat A universe; we estimate that 
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Fig. 3. — Likelihood contours in the h — Qm plane, derived from 
measurements of the BAO observed in the SDSS combined with 
constraints from other cosmological data. The intensities are sep- 
arated by — 21n£ = 2.3, 6.0, 9.2, corresponding to two-parameter 
confidence of 68%, 95% and 99% for a Gaussian distribution. The 
blue shaded region shows the constraints for flat A cosmologies 
combining the SDSS BAO data with a low redshift constraint on 
h = 0.72 ± 0.08 for flat A-cosmologies. The green region shows the 
combination of the SDSS constraint with the 3-year WMAP con- 
straints on U° A f 75 h and Q b h 2 , again for flat A-cosmologies. The 
overlaid red contours were calculated by instead combining with a 
constraint on Qm^ 2 from the peak heights in the CMB together 
with the constraint on Q^h 2 . This relaxes the assumption of a flat 
A-cosmology, by removing the dependence on the distance to the 
last scattering surface. There is still a dependence on the comov- 
ing distance-redshift relation over the survey, but observations of 
type la supernovae constrain this to be close to that expected for 
a flat A-cosmology. The degeneracies in the h — Qm plane induced 
by different co s molog ical observations are discussed at length in 
ITegmark etaf] ||2006T ). 

the allowable residual effect on VLm is at most ±0.02. 

As an alternative to direct h measurement, we con- 
sider the combination of our low-z, redshift-space mea- 
surement of the BAO scale with the angular scale of the 
acoustic oscill ations measured in the CMB. The 3-year 
WMAP data (jSpergel et all 120061) yields the constraint 
Q,°M 75 h — 0.4921° ™? after marginalizing over the scalar 
spectral index and the baryon density. Combined with 
our BAO measurement, this yields Om = 0.2561q'q24; 
and h = 0.7091q'q?,7- For simplicity, we have made the 
good but not perfect approximation that the WMAP 
baryon density constraint, lOOfik/i 2 = 2.233to'ogi; i s m ~ 
dependent of the acoustic scale constraint. Relaxing the 
baryon density constraint by a factor of four makes essen- 
tially no difference to either the best-fit SIm or the error 
bars. This second determination of Qm relies more heav- 
ily on the assumption of a flat A universe, since space 
curvature or a different dark energy component would 
change the distance to the last scattering surface and 
hence the angular scale of the CMB acoustic peaks. 

We can weaken the dependence on the assumption of 
a flat A cosmology by instead combining our low redshift 
BAO measurement with the WMAP constraint on the 
matter density, £1a/^ 2 = 0.12681q;qo95- This constraint 
comes from the relative heights of the CMB acoustic 
peaks, so the physics that underlies it is somewhat dif- 
ferent, and somewhat more degenerate with the parame- 
ters of the inflationary fluctuation spectrum. The depen- 
dence on the assumption of a flat A cosmology is weaker 
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because the relative CMB peak heights only depend on 
the distance to the last scattering surface through the 
weak, and predominantly large-angle, ISW effect. Addi- 
tionally, we have already argued that observations of type 
la supernovae constrain the distance-redshift relation for 
the low-redshift galaxies to be close to a flat A cosmol- 
ogy (leaving a residual effect Afijvf at most ±0.02). The 
direction of the Oa/Zi 2 constraint in the h — Q,m plane is 
more complementary to the low redshift BAO observa- 
tions for measuring Qm compared with the CMB peak 
positions, and this combination yields Qm = 0-256^ 'o22' 
again marginalising over the WMAP baryon density con- 
straint. 

5. DISCUSSION 



In a separate paper, iTegmark et al.1 ((2006) carry out 
full multi-parameter fits to the LRG power spectrum and 
WMAP data over a broader CDM model space, draw- 
ing on both the BAO measurement and the broadband 
shape of the power spectrum. The more focused analysis 
presented in this paper is complementary, obtaining con- 
straints with minimal dependence on detailed cosmologi- 
cal assumptions. Within the context of flat A models, the 
measurement of the matter density, Om = 0.256^0.024' 
from the locations of the acoustic oscillations (the sec- 
ond combination above) is especially "clean," relying on 
a single piece of well understood physics. In particular, 
we have decoupled constraints from the BAO with con- 
straints from the overall shape of the p ower spectrum 
on the same scales. As demonstrated in iPercival et al.1 
(2006), the overall shape of the power spectrum, even 
on scales where the matter clustering has not deviated 
strongly from a linear CDM model is dependent on the 
galaxy population studied. It is therefore imperative to 
accurately model galaxy bias before robust cosmological 
constraints can be derived from such observations. By 
simply considering the BAO in this paper, we avoid this 
complexity. 

Given the different physics probed, the consistency 
of the three Om estimates calculated with different 



additional information, itself provides support for the 
assumptions of a flat universe with a cosmological 
constant. With larger samples and a wider redshift 
range, the BAO "standard ruler" can be used to test 
these assumptions at high pr ecision through the comov- 
ing distance-redshift re lation ([Bond fc Glazebrookll2003t 
ISeo fc Eisensteinl 12003). Forthcoming surveys have been 
designed exploit this effect. It is clear that we are en- 
tering an era where BAO offer an extremely attractive 
route to strong, robust cosmological constraints. 
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